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Abstract. The sulphur chemistry in nine regions in the earliest stages of high-mass star formation is studied 
through single-dish submillimeter spectroscopy. The line profiles indicate that 10-50% of the SO and SO2 emission 
arises in high-velocity gas, either infalling or outflowing. For the low-velocity gas, excitation temperatures are 
25 K for H2S, 50 K for SO, H2CS, NS and HCS+, and 100 K for OCS and SO2, indicating that most observed 
emission traces the outer parts (T < 100 K) of the molecular envelopes, except high-excitation OCS and SO2 
lines. Abundances in the outer envelopes, calculated with a Monte Carlo program, using the physical structures 
of the sources derived from previous submillimeter continuum and CS line data, are ~10-* for OCS, ~10-^ for 
H2S, H2CS, SO and SO2, and ~10"^° for HCS+ and NS. In the inner envelopes (T > 100 K) of six sources, the 
SO2 abundance is enhanced by a factor of ~100-1000. This region of hot, abundant SO2 has been seen before 
in infrared absorption, and must be small, <0'.'2 (180 AU radius). The derived abundance profiles are consistent 
with models of envelope chemistry which invoke ice evaporation at T ~ 100 K. A major sulphur carrier in the 
ices is probably OCS, not H2S as most models assume. Shock chemistry is unlikely to contribute. The source-to- 
source abundance variations of most molecules by factors of ~10 do not correlate with previous systematic tracers 
of envelope heating. Without observations of H2S and SO lines probing warm ( 100 K) gas, sulphur-bearing 
molecules cannot be used as evolutionary tracers during star formation. 

Key words. ISM: molecules - Molecular processes - Stars: Circumstellar matter; Stars: formation 



1. Introduction 

Spectral line surveys at submillimeter wavelengths 
have revealed considerabl e chemical differences between 
star-f or ming regions ( e. g..lBlake et al.lll987HSutton et al. 
1995| ISchilke et all llQQjt iHelmich fc van Dishoeck 
1997t iHatchell et al.l Il998at iNummelin et all I2OO0I: "see 



van 



DishoeclJ hmf for a complete list). While these 



differences indicate activity, the dependence of molec- 
ular abundances on evolutionary state and physical 
parameters is poorly understood. Better insight into this 
relation would be valuable for probing the earliest, deeply 
embedded phases of star formation where diagnostics at 
optical and near-infrared wavelengths are unavailable. 
This is especially true for the formation of high-mass 
stars, for which the order in which phenomena occur is 
much less well understood than in the low-mass case, and 
for which the embedded phase is a significant fraction 
(«10%) of the total lifetime of ^10^ yr. 

Sulphur-bearing molecules are attrac tive as candi- 
date tracers of early protostel lar evolution l)Hatchell et alJ 
Il998bt iBuckle fc Fulled l2n0,1^ . In models of "hot cores" 
^ICharnlevlll997^ . where the chemistry is driven by the 
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evaporation of icy grain mantles due to protostellar heat- 
ing, the abundances of sulphur-bearing molecules exhibit 
a characteristic variation with time. After its release from 
the grain mantles, H2S is converted into SO and SO2 
in -10^ yr, and further into CS, H2CS and OCS af- 
ter ~10^ yr. This type of model has had some suc- 
cess in reproducing the chemistry of star-forming regions 
iLangcr ct al. 2000). 

Besides thermal heating, shocks are potentially rel- 
evant to the chemistry in star-forming regions. Young 
stars drive outflows which expose their surroundings to 
the effects of shocks. The submillimeter lines of SO2 
in Or ion are very broad, and likely arise in the out- 
flow ijSchilke et all Il997^. Models of chemistry behind 
interstellar shocks jMitchel]! Il984l: iLeen fc Gra^ Il988l: 
IFmeau des Forets et alJ Il993|) predict enhancements of 
H2S, SO and SO2 on time scales of ~10^ yr. Thus, both 
by ice evaporation and by shock interaction, sulphur could 
act as a clock on the time scale of 10^ yr relevant for the 
embedded phase of star formation. 

This paper uses submillimeter observations of sulphur- 
bearing molecules to test models of chemistry during star 
formation, and attempts to order these regions chronolog- 
ically. The sources are nine regions of high-mass star for- 
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Fig. 1. Two examples of spectral settings toward eight of our sources, showing differences in the relative and in the 
absolute strengths of lines of sulphur-bearing molecules. The top and bottom frequency scales are the two receiver 
sidebands. The data for NGC 6334 IRSl have been divided by 10 and 3 respectively, for clarity. The non-labeled lines 
are due to non-sulphur-bearing molecules, of which the strongest are labeled in the AFGL 2591 panels. 



mation with luminosities 1 x 10^ — 2 x 10^ Lq at distances 
of 1 — 4 kpc (Table Originally selected for their mid- 
infrared brightness, the sources were found to have large 
envelope masses (40-1100 Mq within r = 0.15 — 0.36 pc) 
based on submillimeter data I van der Tak et all bOO Ob) . 
Together with their weak radio emission and their strong 
molecular outflows, these masses indicate very early evo- 
lutionary stages, probably preceding the formation of hot 
cores. Two of our sources may be more evolved than the 
others: NGC 6334 IRSl, which has a very rich, 'hot core'- 
type submillimeter spectrum l|McCutcheon et alJ Eooot 
Thorwirth et al., in prep.), and NGC 7538 IRSl, which 
has strong free-free emission. 

Obse rvations of CH.sOH lines towards some of these 
sources (jvan der Tak et a'l]l2000a|) indicate an increase in 
the CII3OH abundance by a factor of 100 in the warm 
inner envelope, which is likely due to evaporation of 
CHaOH-rich ices. The excitation of CII3OII, as well as 
that of CO and C2H2, correlates well with other tempera- 
ture tracers, such as the gas/solid ratios of CO2 and II2O, 



the 45/ 100 /im colour and the fraction of heat e d solid 
"CO9 llBoogert et al.l l2000t Ivan der Tak et alJ l2000bt 
iBoonman et alJl2003l) . Since these quantities also corre- 



late with the ratio of envelope mass to stellar mass, their 
variation likely reflects evolution through the dispersal of 
the protostellar envelopes. These findings make the sam- 
ple a good testbed for theories of chemical evolution dur- 
ing the earliest deeply embedded phases of high-mass star 
formation. 



2. Observations 

The data presented in this paper are part of a targeted 
spectral line survey of embedded massive stars, aimed 
to track chemical evolution during the earliest stages 
of high-mass star formation. Fifteen frequency settings 
have been observed, which also contain many lines of sul- 
phurless i nolecules. Some of these data hav e been pub- 
lished by Ivan der Tak et all l|l999t ^Q^Q^. These ob- 
servations were performed with the 15-m James Clerk 
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Maxwell Telescope (JCMT)^ on Mauna Kea, Hawaii be- 
tween 1995 and 1998. The beam size (FWHM) and main 
beam efficiency of the antenna were 18" and 64 — 69% at 
230 GHz and 14" and 58-63% at 345 GHz. The frontends 
were the receivers A2, B3 and B3i; the backend was the 
Digital Autocorrelation Spectrometer, covering 500 MHz 
instantaneous bandwidth. Pointing was checked every 2 
hours during the observing and was always found to be 
within 5". To subtract the atmospheric and instrumental 
background, the chopping secondary was used with 180" 
offsets. Total integration times are 30-40 minutes for each 
frequency setting. Figure shows examples of the data 
towards all but o ne of our sources; the data on W 3 IRS5 
were presented bv lHelmich fc van Dishoecd lll997l) . 

Additional observations of H2S and OCS in the 130- 
180 GHz window were carried out with the 30-m tele- 
scope of the Institut de Radio Astronomic Millimetrique 
(IRAM)2 on Pico Veleta, Spain, in August 2002. The 
frontend was the facility receiver C150 and the backend 
the Versatile Spectral Assembly (VESPA) autocorrela- 
tor. Integration times were 10-20 minutes using frequency 
switching with a throw of 7.3 MHz. Data were calibrated 
onto Tmb scale by multiplying by 1.43, the ratio of for- 
ward and main beam efficiencies. The beam size is 15" at 
these frequencies. 

The H2S 393 GHz line was observed in May 1995 at the 
10.4-m Caltech Submillimeter Observatory (CSO)^, with 
the facility receiver as frontend and the 50 MHz AOS as 
backend. The telescope has a beam size of 25" and a main 
beam efficiency of 60% at this frequency. 

Reduction was carried out using the CLASS package 
developed at IRAM. Linear baselines were subtracted and 
the spectra were smoothed once and calibrated onto Tmb 
scale. Only the frequency switched 30m data required 
higher order polynomials to fit the baseline. The final spec- 
tra have a resolution of 0.3 — 1.5 km s^^ and rms noise 
levels (in Tmb) of 20 - 30 niK for the 150 and 230 GHz 
bands, and 30 — 50 mK for the 345 GHz band. Multiplying 
these values by the line width (Table yields the uncer- 
tainties of the line fluxes in Table|31 Although the absolute 
calibration is only correct to «30%, the relative strength 
of lines within one frequency setting is much more accu- 
rate. 

Li ne identification wa s performed using the JPL cat- 
alog ijPickett et al.lll998l) '^. A matching frequency is the 
main criterion for assignment, but the upper energy level 
of the transition and the complexity of the molecule are 



Table 2. Observed transitions 



^ The James Clerk Maxwell Telescope is operated by the 
Joint Astronomy Centre, on behalf of the Particle Physics 
and Astronomy Research Council of the United Kingdom, 
the Netherlands Organization for Scientific Research and the 
National Research Council of Canada. 

^ IRAM is an international institute for research in mil- 
limeter astronomy, cofunded by the Centre National de la 
Recherche Scientifique (France), the Max Planck Gesellschaft 
(Germany), and the Instituto Geografico Nacional (Spain). 

3 The CSO is supported by NSF grant AST 99-80846. 

* http://spec.jpl.nasa.gov 



Molecule 


Transition 




Frequency 
(MHz) 


(K) 


Telescope Bf 


H2CS 




7l6 ^ 


615 


244047.8 


60.0 


JCMT 


H2CS 


JKpK„ ~ 


IO19 ^ 9l8 


348531.9 


105.2 


JCMT 


HCS+ 


J = 6 


5 




256027.8 


43.0 


JCMT 


HCS+ 


J = 8 ^ 


7 




341350.8 


73.7 


JCMT 


H2S 


JKpKa ~ 


lio ^ 


loi 


168762.8 


27.9 


IRAM 


H2S 


JXpKa ~ 


2ii -» 


2o2 


393450.5 


73.6 


CSO 


H2S 


JXpRa ~ 


220 ~* 


2ii 


216710.4 


84.0 


JCMT 


NS 


'ni/2 J ■- 


= 11/2 


-+ 9/2 


253969 " 


39.9 


JCMT 


NS 


'ni/2 J: 


= 15/2 


13/2 


346220 


71.0 


JCMT 


OCS 


J = 13 ^ 12 




158107.4 


53.1 


IRAM 


OCS 


J = 28 ^ 27 




340449.2 


237.0 


JCMT 


OCS 


J = 30 ^ 29 




364749.0 


271.4 


JCMT 


oc^^s 


J = 18 - 


* 17 




213553.1 


97.4 


JCMT 


oc^^s 


J = 22 -» 21 




260991.8 


144.1 


JCMT 


so 


Nj = 56 






219949.4 


35.0 


JCMT 


so 


Nj = 66 


^ 55 




258255.8 


56.5 


JCMT 


so 


Nj = 67 


^ 56 




261843.7 


47.6 


JCMT 


so 


Nj = 87 


^ 76 




340714.2 


81.2 


JCMT 


so 


Nj = 88 


77 




344310.6 


87.5 


JCMT 


so 


Nj = 89 


^ 78 




346528.5 


78.8 


JCMT 


3*S0 


Nj = 88 


^ 77 




337582.2 


86.1 


JCMT 


S02 


JXpKa ~ 


524 


4l3 


241615.8 


23.6 


JCMT 


S02 


JXpKa ~ 


533 -> 


422 


351257.2 


35.9 


JCMT 


S02 


JXpKa ~ 


533 -> 


524 


256246.9 


35.9 


JCMT 


S02 


JXpKa = 


551 — > 


642 


345149.0 


75.1 


JCMT 


S02 


J Kp Ko — 


'-'33 ' 


624 


254280.5 


41.4 


JCMT 


S02 


JXpKa ~ 


735 


726 


257100.0 


47.8 


JCMT 


S02 


JXpKa ~ 


835 — > 


826 


251210.6 


55.2 


JCMT 


S02 


JXpKa ~ 


937 


928 


258942.2 


63.5 


JCMT 


S02 


JXpKa = 


1064 - 


^ II57 


350862.7 


138.8 


JCMT 


S02 


JXpKa ~ 


111. 11 


lOo.io 


221965.2 


60.4 


JCMT 


S02 


JXpKa ~ 


II39 - 


112,10 


262256.9 


82.8 


JCMT 


S02 


JXpKa ~ 


13l,13 


— > 12o,12 


251199.7 


82.2 


JCMT 


S02 


JXpKa ~ 


14o,14 


^ 13l,13 


244254.2 


93.9 


JCMT 


S02 


JXpKa ~ 


156,10 


-» 165,11 


253956.6 


198.6 


JCMT 


S02 


JXpKo = 


164,12 


163,13 


346523.9 


164.5 


JCMT 


S02 


JXpKa ~ 


18l,l7 


18o,18 


240942.8 


163.1 


JCMT 


S02 


JXpKa ~ 


184,14 


^ 183,15 


338306.0 


196.8 


JCMT 


S02 


JXpKa ~ 


19l,19 


— > 18o,18 


346652.2 


168.1 


JCMT 


S02 


JXpKo ~ 


20i,i9 


192,18 


338611.8 


198.9 


JCMT 


S02 


JXpKo ~ 


222,20 


22i,2i 


216643.3 


248.4 


JCMT 


S02 


JXpKa = 


242,22 


233,21 


348387.8 


292.7 


JCMT 


S02 


JXpKa ~ 


253,23 


^ 252,24 


359151.2 


320.9 


JCMT 


S02 


JXpKa ~ 


282,26 


— > 28l,27 


340316.4 


391.8 


JCMT 


3*S02 


JXpKa ~ 


744 — > 


735 


345519.7 


63.6 


JCMT 


3*S02 


JXpKo = 


1046 - 


^ 1037 


344245.4 


88.4 


JCMT 


3*S02 


JXpKa ~ 


II39 - 


^ 112,10 


253936.3 


81.9 


JCMT 


3*S02 


JXpKa ~ 


152,14 


14l,13 


358988.0 


118.7 


JCMT 


3*S02 


JXpKa ~ 


19l,19 


— > 18o,18 


344581.1 


167.4 


JCMT 


3*S02 


JXpKa ~ 


299,21 


— > 308,22 


257466.7 


590.7 


JCMT 


" Blend of the F = 


7-^6 


, F = Q 


5 and F 


= 5 ^ 


4 



hyperfine components 
''Blend ofthef = 9^ 
components 



F = 8 ^ 7 and F = 7 ^ 6 hyperfine 
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Table 1. Source sample.'^ 



Source 


R.A. (1950) 


Dec. (1950) 


L 


d 


iV(H2) 




(h m s) 


(° ' ") 


(lO"* Lo) 


(kpc) 


(lO^^* cm-2) 


W3 IRS 5 


02 21 53.1 


+61 52 20 


17 


2.2 


2.3 


MonR2 IRS3 


06 05 21.5 


-06 22 26 


1.3 


0.95 


1.8 


NGC 6334 IRSl 


17 17 32.0 


-35 44 05 


11 


1.7 


11.3 


W33A 


18 11 43.7 


-17 53 02 


10 


4 


6.2 


AFGL 2136 


18 19 36.6 


-13 31 40 


7 


2 


1.2 


AFGL 2591 


20 27 35.8 


+40 01 14 


2 


1 


2.3 


S140 IRSl 


22 17 41.1 


+63 03 42 


2 


0.9 


1.4 


NGC 7538 IRSl 


23 11 36.7 


+61 11 51 


13 


2.8 


6.5 


NGC 7538 IRS9 


23 11 52.8 


+61 10 59 


4 


2.8 


3.3 



° Luminosity an d dista nce from iHenning et al.1 lll992ll and IGiannakopoulou et alJ (Il997ll (MonR2 IRS3) and from 
Ivan der Tak et all l)2000lJ) (other sources); H2 column density from this work (MonR2 IRS3) and from Ivan der Tak et all 
(t2000b1 (other sources). 



also considered. Table [21 has details about the observed 
transitions. 

For each detected feature, a satisfying match could 
be found in the catalog within 0.5 MHz, except for the 
H2CS IO19 ^ 9i8 line at 348.5 GHz, whose residual is 
2 MHz. Labo ratory spectroscop y of H2CS only exists up 
to 250 GHz ((Beers et al.lll972l) . and the analysis in the 
JPL catalog only includes pure rotation and centrifugal 
distortion terms. At higher frequencies and for higher J- 
values, higher order terms become important, which may 
be why the catalog prediction is inaccurate for this line. 
Our data indicate a frequency of 348534.2 MHz. 

3. Results 

3.1. Line profiles 

We have fitted the detected lines with Gaussian pro- 
files. Table reports the observed fluxes and Table 01 
the widths of the lines. While most lines could be satis- 
factorily fitted with a single Gaussian, the lines of SO, 
which are the strongest lines that we have observed, 
required two Gaussians for a good fit to their profiles 
(Fig. [2 Table 0)). One component agrees in position and 
width with the va l ues me asured for C^^O and C'^^S by 
Ivan der Tak et all l|2000b|) and is attributed to a static 
envelope. Table [SI reports similar measurements for the 
source MonR2 IRS3 which van der Tak et al. did not 
study. The other component has 2.8 times as large a 
width on average. This high- velocity emission may trace 
outfiowing or infalling motions. The contributions from 
high- and low- velocity gas to the SO line fiuxes are about 
equal, unlike CS, where the env elope carries 90% of the 
fluxes l|van der Tak et alJl2000bl Table [SJ. The infrared 
absorption and submillimeter emission lines of CO to- 
wards these sources also show multiple velocity compo- 
nents jMitchefl et all FlQQlt iMitchell k Haseeawai Il99lt 
{Mitchell et al.ll992|) . However, the velocities do not match 
in detail, as expected because of differences in spectral 
and angular resolution between infrared and submillime- 



ter data on one hand, and excitation differences between 
CO and SO on the other. 

The fit of two Gaussians to the SO lines is excel- 
lent, which is not always the case for CO. An alternative 
method is to use area integrals, which do not attribute 
any low- velocity gas to the high- velocity component. Our 
approach may overestimate the high- velocity contribution 
to the line flux by factors up to 2. 

The ratios of the high- to low- velocity contributions to 
the SO lines are ^ 1 for the lines in the 230 GHz band, 
and Z 1 for the lines in the 345 GHz band. This differ- 
ence could be due to the higher temperatures and den- 
sities required for excitation of the 345 GHz lines, or to 
the smaller beams in which they are observed. Since the 
high-velocity wings on the CS lines in these sources also 
become more pronou nced in smaller be ams, independent 
of transition fvan der Tak et allbOOOhl) . the broad com- 
ponent probably is more compact than the envelope, but 
not necessarily warmer or denser. The high- velocity gas in 
these sources thus seems to be confined to small radii. This 
situation differs from that in molecular clouds, where line 
width increases with size, and from that in many molec- 
ular outfl ows, where velo cities usually increase with ra- 
dius (e.g.. lLee et aLlEoO^) . However, while CO data show 
arcminute-sized outflows in our sources, higher velocities 
naturally occur at smaller radii in the case of infall. For 
example, gas at 1,000 AU (1" at 1 kpc) in free fall onto a 
10 M0 object would have Vinf = 4 km s^^, consistent with 
the observed width of the broad component (Table ^J. 

The broad component peaks at slightly more 
blueshifted velocities (by 0.17 km s~^ on average) than 
the envelope, a situation found before for CS and H2CO 
l)van der Tak et al.ll200 0bV In the case of outflow, obscu- 
ration by Ay = 10"* mag of dust could explain the ab- 
sence of redshifted wings. Since some redshifted SO is seen 
(Fig. [2]), the requirement relaxes to several 1000 mag. Of 
the sources discussed here, only NGC 6334 IRSl may have 
such a large amount of dust in its envelope (Table the 
other sources may need circumstellar disks seen face-on. In 
the case of infall, a mixture of absorption and emission is 
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Table 3. Observed line fluxes / TuBdV (K km s^^) from Gaussian fits, and Icr upper limits. Dots denote missing 
data. 



Transition 


W3° 


AFGL 
2591 


NGC 7538 
IRSl IR89 


MonR2 
1R83 


8140 
1R81 


NGCd334 
IKSl 


AFGL 
2136 


W33A 


H2CS 7l6 — > 615 


0.8 


1.4 


3.8 


1.6 


0.9 


1.5 


10.6 


0.3 


2.9 


H2CS IO19 ^ 9l8 


<0.3 


1.1 


3.5 


1.2 


0.3 


0.8 




1.0 


3.6 


HCS+ 6^5 


0.3 


















HCS+ 8^7 


0.3 


0.6 


1.6 


0.8 




0.7 


16.6 


0.3 


1.1 


H2S lio — > loi 


9.4 


13.1 


15.4 


7.7 


7.6 


14.3 




6.7 


17.8 


H2S 2ii — > 2o2 


2.6 




1.6 














H2S 220 — > 2ii 


2.2 


0.8 


1.9 


0.4 


<0.2 


0.6 


22.8 


<0.1 


1.9 


NS 11/2 ^ 9/2 




0.9 


0.9 


0.5 


<0.4 


0.5 


15.9 


<0.2 


2.2 


NS 15/2 13/2 




<0.8 


<1.2 








11.6 




2.8 


DCS 13 ^ 12 




0.4 


0.8 


0.3 




<0.3 








OCS 28 ^ 27 


0.9 


0.3 














2.6 


OCS 30 ^ 29 


<0.4 


0.3 
















OC^*S 18 ^ 17 














4.8 




0.4 


OC^*S 22 ^ 21 


<0.1 












11.6 






SO 56 ^ 45 


17.4 


9.9 


8.4 


1.3 


3.7 


10.8 




2.0 


3.0 


broad 


11.2 


4.5 


6.5 


5.2 


3.6 


3.7 




1.2 


5.1 


SO 66 ^ 55 


26.1 


5.3 


5.8 


2.5 


3.5 


4.0 


26.1'= 


0.4 


3.1" 


broad 


22.1 


2.5 


5.7 


1.5 


1.4 


2.6 




0.6 




SO 67 56 


14.7 


7.3 


10.0 


3.0 


4.8 


5.7 




1.5 


3.5 


broad 


13.2 


3.7 


7.3 


3.0 


3.1 


4.7 




1.3 


5.6 


SO 87 ^ 76 


17.9 


3.7 


6.5'' 


1.3 




2.6 


14.5 


<0.2 


1.8^= 


broad 


26.8 


5.1 




2.4 




5.5 


23.4 


<0.2 




SO 8g ^ 77 


31.7 


3.4 


9.5 


4.7'' 




2.6 


47.3^= 


1.5 


4.1 


broad 


28.8 


4.5 


6.7 






5.0 




0.5 


6.4 


SO 89 ^ 78 


27.9 


5.6 


8.3 


1.5 


3.6 


3.9 


26.4 


1.4^= 


2.2 


broad 


58.1 


5.1 


5.9 


3.9 


5.7 


6.0 


34.4 




6.7 


^*SO 88 77 


18.9 


2.7 


0.9 


<0.3 




<0.4 


19.6 




<0.7 


SO2 524 4i3 


8.2 


2.3 


2.0 


0.7 


1.0 


1.1 


6.8^* 






SO2 533 422 


12.2 
















2.9 


SO2 533 ^ 524 


6.9 










<1 








SO2 551 642 


<0.3 












7.1 






SO2 633 624 


20.8 


3.0 


1.7 


0.6 


1.6 


1.1 


12.1 


<0.1 


2.2 


SO2 735 ^ 726 


10.2 


1.7 


1.1 


<0.4 


0.7 


0.5 






1.0 


SO2 835 ^ 826 


12.7 


1.7 


0.9 


0.4 


0.7 


0.7 


5.9 


<0.2 


1.4 


SO2 937 — > 928 


9.6 


2.4 


1.5 


<0.2 




<0.5 






3.7 


SO2 1064 ^ II57 


3.4 


1.7 


<0.4 












0.8 


SO2 111, 11 10o,io 
















0.3 




SO2 II39 112,10 


















1.6 


SO2 13l,l3 ^ 12o,12 


19.0 


2.7 


2.2 


0.7 


1.5 


1.5 


8.8 


<0.2 


1.7 


SO2 14o,14 ^ 13l,13 


12.2 


2.8 


1.7 




1.3 


1.3 




0.3 


1.9 


SO2 156,10 165,11 




0.6 
















S02 164,12 — > 163,13 


29.5 














0.6 


2.0 


SO2 18l,l7 18o,18 


7.7 


1.0 














0.7 


S02 184,14 — *■ 183,15 


13.3 


2.9 






1.2 








1.3 


SO2 19l,l9 18o,18 


31.7 


2.6 


<0.4 


<0.4 




<0.7 


7.0 


0.4 




SO2 20l,l9 ^ 192,18 


16.0 


2.7 
















SO2 222,20 22i,2i 


5.5 


0.8 










1.7 




0.5 


SO2 242,22 233,21 


6.1 


3.2 


<0.6 














SO2 253,23 ^ 252,24 


4.1 


2.8 


1.4 












1.0 


SO2 282,26 — * 28l,27 


11.3 
















0.6 


^*S02 744 735 


2.4 


1.0 
















3*S02 1046 ^ IO37 


2.5 


1.1 
















^*S02 II39 — > 112,10 




1.1 














0.8 


^''S02 152,14 -> 14l.l3 


0.7 


1.1 














0.8 


^''S02 19l,l9 18o,18 


5.8 


1.2 


<0.4 


<0.4 




<0.4 








='''802 299,21 ^ 308.22 


<0.2 
















0.3 



" From lHelmich fc van Dishoeckl {li^, except SO 
Attributed to envelope based on center velocity 
Attributed to envelope based on line width 
Difficult fit due to partial blend with other lines 
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J (km s ' 



60 -50 -40 -30 -30 
^LSR (km s"') 



Table 5. Observations of CS, 
wards MonR2 IRS3. 



C^-^S and C"0 lines to- 



Line 


Frequency 
(MHz) 


Tmb 
(K) 


AV 

(km s ) 


J TMBdV" 
(K km s ) 


CS 5 ^ 4 


244935.6 


9.2 


2.7 


27 




broad 


0.2 


8.0 


1.9 


CS 7^ 6 


342883.0 


6.6 


2.5 


17.8 




broad 


0.3 


8.0 


2.4 


CS 10 ^ 9 


489751.0 


2.0 


2.7 


5.6 


C^*S 5 ~* 4 


241016.2 


0.6 


2.3 


1.4 


C^*S 7^6 


337396.7 


0.4 


2.6 


1.1 


C^'^O 3^2 


337060.9 


2.3 


2.6 


6.4 



Fig. 2. Observed velocity profiles of the SO 65-55 and 87- 
Te lines in the sources AFGL 2591 and W3 IRS5, showing 
the low- and high-velocity components. The dashed lines 
are our fits to the individual components and the dotted 
line is the sum of the two components. 



expected at redshifted velocities l|Choill200^ . which may 
partly cancel each other if unresolved, to give the impres- 
sion of predominantly blueshifted emission. We conclude 
that both infall and outflow are plausible explanations for 
the high- velocity gas, and recommend interferometric ob- 
servations of SO to decide between these options. 

For most sources in our sample, the width of the SO2 
lines is between the values for the low- and high-velocity 
components. The only source where the SO2 lines are 
strong enough to fit two components is W3 IRS5, where 
the line profiles do not show wings. However, most o f these 
data, taken from iHelmich &: van Dishoeckl \\m% . were 
observed with an AOS as backend, and have lower spec- 
tral resolution than our data. It seems therefore plausible 
that, like for SO, about half of the SO2 emission arises in 
high- velocity gas. 

For the other molecules, the signal to noise ratio is 
not high enough to disentangle multiple velocity compo- 
nents. This paper assumes that these lines originate in 
low-velocity gas, which given our results for SO may over- 
estimate column densities and abundances (SS 13.31 1^ by 
factors of up to 2. 



ture. Another complication is that the lines are measured 
in beams of unequal size, so that somewhat different vol- 
umes of gas are probed. Despite these caveats, the rotation 
diagram provides a useful first estimate of excitation con- 
ditions and molecular column densities, and a stepping 
stone towards a more sophisticated analysis (§^. 



Table 6. Rotation temperatures (K). Numbers in brackets 
denote uncertainties. 



Source 


SO 




SO2 


CH3OH" 




Narrow 


Broad 






W3 IRS5 


>100 


>100 


154(14)" 


73(32) 


W33A 


54(18) 


96(64) 


124(12) 


155(19) 


AFGL 2136 


59(24) 


34(9) 


>200 


143(12) 


AFGL 2591 


31(6) 


>100 


185(28)'' 


163(9) 


S140 IRSl 


26(4) 


>100 


56(16) 


41(4) 


NGC 7538 IRSl 


65(27) 


61(25) 


47(12) 


189(30) 


NGC 7538 IRS9 


48(17) 


45(15) 




29(2) 


MonR2 IRS3 


47(20) 


>100 


125(35) 


110(20) 


NGC 6334 IRSl 


75(72) 




45(15) 


213(37) 



" Better fit with two components: T^ax= 70(14) & 254(20) K. 
'' Better fit with two components: riot= 47(11) fc 197(67) K. 

Fro m JCMT data llvan der Tak et al-2000a : .Boonman et alJ 
I1999I) . 



3.2. Excitation temperatures 

For SO and SO2, the number of detected hues is large 
enough to construc t rotation diagrarn s. Th is method, de- 



scribe d in detail bv lBlake et al.l 1 19871) and'Hel mich et all 



l|1994|) . assumes that the lines are optically thin and that 
the molecular excitation can be described by a single 
temperature, the rotation temperature. If radiative decay 
competes with collisional excitation, the rotation temper- 
ature lies below the kinetic temperature. If the lines arc 
optically thick, the rotation diagram underestimates the 
column density and overestimates the excitation tempera- 



TableElpresents the results. In the case of SO, separate 
fits were made to the low- and high-velocity components. 
For NGC 7538 IRS9, the detected SO2 Hues do not cover a 
large enough range of energy levels to constrain T^^t ■ The 
derived temperatures are 30. ..70 K for SO and 50. ..190 K 
for SO2, reflecting the available range of energy levels. 
Lower limits to Tiot are poor fits to the data, probably 
caused by optically thick lines. Indeed, the SO/'^'^SO and 
S02/'^^S02 line ratios of 1.3... 10 are below the isotopic 
abundance ratio of ^^S/^'*S=:22, which indicates optical 
depths of a few in the SO and SO2 lines. 
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Table 4. Widths AV (km s^^) of the emission hne components, averaged over transition. Numbers in brackets denote 
uncertainties in units of the last decimal, except when only one transition was observed. Dots denote missing data. 
For SO, values for the narrow and the broad components are reported, whose positions are also given. 



Source 


CS 




SO (narrow) 


SO (broad) 


SO2 


ocs 


H2S 


H2CS 


IICS+ 




Vlsr 


AV 




AV 




AV 












W3 IRS5 


-38.4(3) 


2.7(11) 


-38.7(5) 


4.8(3) 


-39.0(10) 


13.5(14) 


6.6(12) 


3.5(20) 


5.6(6) 


2.9 


1.1(4) 


W33A 


+37.5(9) 


5.4(5) 


+37.8(3) 


3.6(2) 


+36.9(10) 


11.6(11) 


5.0(3) 




2.5 


2.2(7) 




AFGL 2136 


+22.8(1) 


3.1(4) 


+22.9(2) 


2.5(8) 


+23.6(16) 


6.5(13) 


5.0(7) 


4.1(2)" 


5.4 


6.0 


5.4 


AFGL 2591 


-5.5(2) 


3.3(6) 


-6.0(3) 


3.7(2) 


-5.6(9) 


8.6(15) 


5.0(30) 


4.50 


4.1 


5.1(12) 


7.1 


S140 IRSl 


-7.0(2) 


3.3(3) 


-7.0(1) 


2.5(2) 


-8.1(2) 


9.6(22) 


1.9(10) 




3.6 


2.9(8) 


1.0 


NGC 7538 IRSl 


-57.4(5) 


4.1(14) 


-57.5(4) 


3.8(4) 


-58.4(6) 


11.5(39) 


5.1(9) 


3.8(1) 


3.8(2) 


3.2(4) 


3.2 


NGC 7538 IRS9 


-57.2(3) 


4.1(5) 


-57.4(3) 


3.2(8) 


-55.9(5) 


7.8(21) 


3.5(8) 




2.9(1) 


2.9 


2.6 


MonR2 IRS3 


+10.2(1) 


2.4(2) 


+10.1(1) 


3.0(1) 


+10.0(16) 


8.2(4) 


7.0(9) 


3.1 


4.0(5) 


4.4 


3.6 


NGC 6334 IRSl 


-7.4(2) 


5.3(3) 


-7.2(11) 


4.1(1) 


-7.4(13) 


11.1(1) 


5.2(10) 


2.1 


4.1(3) 


3.5(1) 


3.5 



" Value for OC^^-'S 



i4 



CO 

O 
C/3 



200 



150 



° 100 



50 



I — <> — I 



50 100 150 200 
r,„t(CH30H) (K) 

Fig. 3. Rotation temperatures of SO (dots) and SO2 (tri- 
angles) versus values for CH3OH. 



The last column of Table El lists rotation tempera- 
tures of CH3OH, measured previously. Strictly, these num- 
bers are only lower limits to the kinetic temperatures 
in a 15" beam, but they show good correlation with 
Tr^(C-;i ll-?,) which does directly t race the kinetic temper- 
ature ijvan der Tak et a'i]l2000a|) . The rotation tempera- 
tures of SO and SO2 are not clearly correlated with Trot 
(CH3OH) within the uncertainties (Fig. Therefore we 
cannot conclude at this point whether these molecules 
trace the same gas within the envelopes. In any case, dif- 
ferences in beam sizes and optical depth effects in the lines 



make it hard to draw firm conclusions before a more de- 
tailed radiative transfer analysis is presented in § 0] 

For the other molecules, for which only few lines are 
available, the line ratios have been used to constrain Tex- 
Source-averaged line ratios indicate Tcx~50 K for H2CS, 
NS and HCS+, Tex«25 K for H2S and Tex«100 K for OCS. 
These values are only rough estimates, especially if based 
on few sources and/or lines. 



3.3. Column densities 



Table presents beam-averaged molecular column den- 
sities, estimated from the observed line fluxes using the 
excitation temperatures found above. The column densi- 
ties for H2CS and H2S are the sums of the ortho and para 
species, assuming an ortho to para ratio of 3. After cali- 
bration, the adopted Tex is the main source of uncertainty 
in the column densities, especially if only one or two lines 
have been observed. For example, decreasing Tqx from 50 
to 25 K increases column density estimates by factors of 
ss2, while increasing Tex from 50 to 75 K decreases column 
density estimates by ^ 20%. In addition, column densities 
may have been underestimated by factors of a few due to 
nonnegligible optical depth. 

To search for trends of column densities with tem- 
perature, Table \7\ lists the sources in order of increasing 
T(CH30II). The only trend seems to be that the warmest 
source has the largest column density of all molecules, ex- 
cept SO and SO2. In order to investigate whether these 
trends in the column density reflect chemical differences 
between sources or are due to optical depth or beam size 
effects, a more detailed radiative transfer analysis is pre- 
sented in § 01 
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Table 7. Column densities (cm^^) in 15-20" beams. 



Source 


SO (10 




SO2 


OCS 


H2S 


H2CS 


NS 


HCS+ 




Narrow 


Broad 


10^3 


10^^ 


10" 


10" 


10" 


10" 


NGC 7538 IRS9 


3 


4 




0.6 


2 


6 


1 


5 


S140 IRSl 


10 


>6 


7 


<0.8 


3 


7 


1 


6 


W3 IRS 5 


>27 


>11 


530" 


2 


7 


5 




3 


MonR2 IRS3 


5 


>5 


15 




3 


4 


<3 




AFGL 2136 


2 


1 


>3 




2 


4 


<2 


3 


W33A 


4 


9 


26 


9 


4 


10 


4 


6 


AFGL 2591 


10 


>7 


52'' 


1 


3 


6 


2 


5 


NGC 7538 IRSl 


12 


9 


11 


2 


5 


15 


2 


11 


NGC 6334 IRSl 


40 




46 


24 


60 


30 


20 


84 



° Better fit with two components: 7V= 
' Better fit with two components: A'^= 



1.6 and 5.5 xlO'* cm 



3.4. Comparison with infrared data 

iKeane et al observed the i>3 band of SO2 around 

7.3 ^m in absorption towards five of the sources stud- 
ied here, and found column densities of SO2 of a few 
10^^ cm~^, a factor of ^^100 higher than the submillime- 
ter data indicate. In contrast, infrared observations of CO 
and dust yield 3-5 times lower column densities than sub- 
millimeter data, as a result of non-spherical geometry on 
scales ^ 1" l)van der Tak et all2000b'l . Correcting for opti- 
cal depth in the SO2 submillimeter lines increases the col- 
umn density estimate, but only by factors of a few fS I3.2l) . 
which does not explain the discrepancy with the infrared 
values. More likely to be important are chemical effects 
on small scales. The sources are centrally concentrated, 
n(H2)oc r~" with a = 1.0 — 1.5, so that absorption data 
are more sensitive to warm, dense gas at small radii, while 
emission data probe more extended, cooler and less dense 
gas. The infrared excitation temperatures of 225-750 K 
also suggest that the SO2 absorption arises in warm gas. 

The infrared estimates of (SO2) imply optically thick 
submillimeter lines of SO2, but their brightness measured 
with the JCMT is much lower than the infrared Tex- The 
implied limit on the source size is 6*5 = 15"(Tox/Tb)~°'^, 
with To the typical Tmb of 0.1 K. Inferred values of 9s 
range from 0'.'17 for AFGL 2591 to 0'.'32 for MonR2 IRS3, 
corresponding to linear radii of 100 — 250 AU. The lim- 
ited spectral resolution of the infrared data prohibits as- 
signment of the absorbers to either velocity component in 
emission, which however does not affect our conclusion. 

4. Radiative transfer analysis 

4.1. Model description 

To estimate molecular abundances we have used 
the Monte Carlo radia ti ve tr ansfer program by 
iHogerheiide &: van der Ta^ l)2000(l ^. Starting from 
radial profiles of the density and kinetic temperature, this 

^ http:/ /talisker. as. arizona.edu/~michiel/ratran.html 



program solves for the molecular excitation as a function 
of radius. Besides collisional excitation, radiation from 
the cosmic microwave background and thermal radiation 
from local dust are taken into account. The result is 
integrated over the line of sight and convolved with the 
appropriate telescope beam. Observed and synthetic line 
fluxes are compared with a statistic to find the best-fit 
abundance. Initially, abundances were assumed to be 
constant with radius. 

Molecular data were mostly taken from the database 
by Schoier et al. (in prep.)^, which paper also de s cribes 
our extension of the SO2 rate coefficients bv lCreenl l)l99.4l 
to higher energy levels. No rate coefficients are available 
for NS, so its excitation was assumed to be thermalized 
at the temperature of each grid point of the model. The 
population of the ^113/2 state of NS, which lies 225 K 
above the ^^\/-? state, was assum ed to be negligible. In 
the case of H^S. lBall et all 1)1999(1 measured the lio — loi 
rate coefficient in collisions with He at T = 1.36 — 35.3 K . 
Our model uses the rate coefficients bv lCreen et al 1 lll993l) 
for H20-IIe, multiplied by 2.6 based on comparison with 
the Ball et al. data, and scal ed for the different reduced 
mass of the H2S-H2 system. iTurner. (.1996) lists reasons 
why II2S rate coefficients are likely to be larger than II2O 
values. 



4.2. Physical structure 

The temperature and density structure s of all our sources 
but M onR2 IRS3 were derived by Ivan der Tak et all 
ll2000bl) . We have followed the same procedure to in- 
fer the structure of MonR2 IRS3, assuming a power- 
law density structure n — no(r/ro)^". The temperature 
profile is deri ved from the o bserved luminosity of 1.3 x 
10"* L0 (Henning et al.lll993) . whi le m is constrained by 
the su bmillimeter photometry of ICiannakopo^ilou et alJ 
1I1997I) . These dust models solve the grain heating and 
cooling self-consistenly, using opacity model 5 from 



http://www.strw.leidenuniv.nl/ ~moldata 
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lOssenkopf fc Hennind l)l994 . These "0H5" opacities are 
the only ones that yield dust masses consistent with 
C^^O measurements of AFGL 2591, where CO deple- 
tion is known to be negligible l)van der Tak et al]ll999|l . 
To constrain a, the CS and C^'*S line spectrum of 
MonR2 IRS3 was modeled with the Monte Carlo pro- 
gram. Be sides our own . JCMT m easurements (T able 13), 
data from lTafalla et alJ l)l997(l and lChoi et all ^00) were 
used. The best fit was obtained for a = 1.25 (Table |H1), 
with a reduced = 2.33 over 11 degrees of freedom, al- 
though a = 1.0 (x^ = 3.1) and a = 1.5 (x^ = 3.8) also 
give acceptable fits. Radial profiles of 350 /im dust emis- 
sion for our sources indic ate larger values of a, by up to 
0.5, than CS line emission l)Mueller et alJ2002(l . If this dif- 
ference holds more generally, it may reflect a reduced CS 
abundance at small radii, although the models of ii l5.2l do 
not predict that. 

Table 8. Power law model n = no{r/ro)~°' for the density 
structure of MonR2 IRS3. 



Parameter 



Value 



Outer radius ro 
Inner radius 
a 

no (H2) 
CO/H2 
Dust mass 



0.3 pc 
0.001 pc 
1.25 
1.21 X W cm" 
8 X 10"^ 
1.05 M0 



4.3. Molecular abundances 

Tabic inihsts the results of the radiative transfer analysis, 
with the same source ordering as in Tabled Abundances, 
if assumed constant, are ~10^^ for OCS, '--^10^^ for H2S, 
H2CS, SO and SO2, and -10^1° for HCS+ and NS. These 
trends are the same as found above for beam-averaged 
column densities (S I3.3|l : the Monte Carlo analysis does 
not change results by factors of more than a few. Except 
for SO, SO2 and OCS, the source-to-source spread in the 
abundances is less than a factor of 10, and not related to 
temperature. These molecules (CS, H2S, H2CS, NS and 
HCS"*") appear to trace the chemically inactive outer en- 
velope (T < 100 K), at least in the observed transitions. 
However, the high OCS abundances occur all in warm 
sources. 

For SO and SO2, the number of detected lines is large 
enough to investigate their abundances as a function of ra- 
dius. Figure 01 shows that the Monte Carlo models fit the 
SO lines uniformly well over the observed range of energy 
levels for all sources except W 33 A: there is no indication 
for changes in the SO abundance between T = 30 and 
90 K. In the case of W 33 A, some of the observed values 
appear to be higher than the models, which may be due 
to an uncertain high-velocity contribution for these lines. 
Observations of higher-excitation SO lines (TV = 10 — > 9 



near 430 GHz or iV = 11 10 near 473 GHz) would be 
valuable to investigate the effect of evaporating ice man- 
tles on the abundance of SO. 

Figure [S] shows the results of the models for SO2. 
Models with SO2/H2 ~ lO'^ fit the lines with K < 100 K 
within factors of a few. However, these models underpro- 
duce the hues with Eu 1 100 K by factors of > 100. This 
discrepancy is best seen for W3 IRS5, W33A and AFGL 
2591, where the most SO2 hues are detected. For these 
sources, we have run models where the SO2 abundance in- 
creases sharply ('jumps') when T > 100 K. These models 
are shown as triangles in Fig. Eland quantified in Table [Till 
The SO2 abundance increases by a factor of 100 in the case 
of W3 IRS5 and W33A, and by a factor of 1000 for AFGL 
2591. These models reproduce all observed SO2 lines to 
factors ~10. The residuals are considerably bigger than 
the expected error bars of factors of 2... 3, mainly due 
to the uncertain high-velocity contribution, which sug- 
gests that more detailed modeling is needed. However, 
these 'jump' models fit the data much better than models 
where the SO2 abundance is constant. The SO2 abun- 
dances derived from the high-excitation lines are similar 
to those based on the infrared data, indicating that the 
same gas is probed. The data for the other sources are 
consistent with similar abundance jumps, but not enough 
high-excitation SO2 lines have been detected to constrain 
such models. However, the data of AFGL 2136, NGC 7538 
IRSl and MonR2 IRS3 allow estimates of SO2 abundances 
in their inner envelopes, which Table 11111 also lists. Except 
NGC 6334 IRSl, the jumps in the SO2 abundance occur 
all in sources with high r(CH30H), confirming that the 
abundance enhancement is related to high temperatures, 
probably through ice evaporation. 

Table 10. Models with SO2 abundance 'jumps'. 



Source 


(S02/H2)out 

10-8 


(S02/H2)in 

10-^ 


W3 IRS5 


10 


10 


MonR2 IRS3 


1 


12'' 


AFGL 2136 


0.5 


0.8' 


W33A 


1 


1 


AFGL 2591 


2 


2 


NGC 7538 IRSl 


1 


146 


IRAS 16293" 


0.5 


1 



° From lSchoier et alJ EoO: 
^ Estimated value. 



The models discussed so far assumed a static envelope 
and aimed to describe the low- velocity gas. To investigate 
whether the high-velocity gas could be due to infalling 
motions, we have modeled the SO emission from W3 IRS5 
under the assumption of gravitational infall. In this model, 
the infall velocity I^nf varies with radius R as y/GM/R, 
where G is the gravitational constant. The luminosity of 
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Table 9. Molecular abundances and abundance ratios: outer envelope 



Source 


SO 


SOa 


OCS 


CS" 


HaS 


HaCS 


NS 


HCS+ 


SO/ 


SOa/ 


SO/ 


HaS/ 


NS/ 




10-9 


10-9 


10-9 


10-9 


10-9 


10-10 


10-" 


10-"' 


SOa 


CS 


CS 


CS 


CS 


NPP 7^'^8 TR'^Q 


1 
1 






1 


Q 
O 


1 


1 

_L 


p; 
o 


9 


u.uo 


n 1 

U. 1 


u.o 


n nni 

U.UUl 


S140 IRSl 


2 


1 


<2 


5 


2 


7 


1 


5 


2 


0.2 


0.4 


0.8 


0.002 


W3 IRS5 


10 


10 


0.5 


5 


3 


3 




0.2 


1 


2 


2 


2 




MonR2 IRS3 


5 


1 




5 


<2 


7 


<1 




5 


0.2 


1 


<0.8 




AFGL 2136 


1 


0.5 




4 


<1 


3 


<0.5 


1 


2 


0.1 


0.25 


<0.5 


<0.001 


W33A 


1 
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W3 IRS5 , if due to a single star, suggests a stellar mass of 
«30 Mq ijVacca et aLlll996l) . so that T^nf ranges from 0.45 
to 10.66 km s"^ in the model. However, the synthetic SO 
line profiles in 14-18" beams have widths of 5-6 km s^^ 
FWHM, much less than the observed value of 13.5 km s^^. 
Since the other sources have lower luminosities (hence in- 
fall speeds), it seems unlikely that the high- velocity SO 
and SO2 emission is solely due to infall. 



In the case of hot cores, both abundance ratios and ab- 
solute abundances are similar to those observed here. The 
chemical similarity of our sources with hot cores prompts 
us to compare our results with chemical models invoking 
ice evaporation. Our assumption is that source-to-source 
differences in abundances are due to age differences, rather 
than to differences in initial physical or chemical condi- 
tions. 



5. Discussion 

5.1. Comparison to other objects 

Tables and ^| also list molecular abundances in other 
regions. Only the case of the low-mass protostar IRAS 
16293 has been analyzed to the same level of detail as done 
here. To reduce uncertainties introduced by methodologi- 
cal differences, the table contains abundances scaled to the 
value of the commonly observed CS molecule. Using ratios 
also facilitates comparison with chemical models, because 
these are independent of the total amount of sulphur in the 
gas phase, which is not well-understood. Ultraviolet obser- 
vations of diffuse clouds toward C, Oph indicate that sul- 
phur is undepleted from i ts solar abundance of ~2xl0-^ 
ijSavaee &: SembacbllT99^ . However, for dense clouds, it 
is unclear what fraction of sulphur is locked up in dust 
grains. 

The ratios observed in our sources differ considerably 
from those in PDRs, dark clouds and shocks, indicating 
that in those regions, other chemical processes dominate. 
The ratios in IRAS 16293 are at the low end of the values 
found here, suggesting that in low-mass objects, depletion 
of molecules onto grains is more important. 



5.2. Comparison with chemical models 

The sulphur chemistry of hot cores, where ice mantles are 
evaporating off dust grains and start a high -temperature 
chemistry, was studied by ICharnlevI (|l99a). The model 
assumes that the dominant form of sulphur in the ices is 
H2S, which after evaporation is destroyed by H and H3O+ 
to form S and HaS"*". Reactions with OH and O2 make 
first SO and then SO2. Some H2S reacts with into CS, 
which reacts with OH to make OCS. The sulphur chem- 
istry strongly depends on temperature through its con- 
nection with the oxygen chemistry. At T> 300 K, the SO2 
abundance remains factors of 10-100 below that of SO be- 
cause O and OH are driven into H2O. Results also depend 
on the presence of O2 in the ice man tles. Observational 
limits on solid and gaseous O^ l|Vandenbussche et al.ll999l: 
iGoldsmith et alJl2000t IPaeani et alJl2003|) do not firmlv 
rule out either assumption, but do make the model with- 
out O2 ice more plausible. However, this model cannot 
directly be compared with our observations, which probe 
a mixture of warm and cold gas. 

The chemis t ry of protostellar envelopes was modeled 
by iDotv et "all l)2002|) . In the inner warm region, a hot 
core chemistry similar to that of iCharnlevi ((1992) is used. 
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whereas in the cold outer parts, a low-temperature chem- 
istry including freeze-out is adopted. The initial conditions 
of the model depend on temperature to mimic the effect 
of ice evaporation. Although Doty et al. adopt the specific 
tempera ture and density st r ucture of AFGL 2591 as mod- 
eled by Ivan der Tak et all l|2000hl) . their results do not 
change by factors of more than a few for our other sources, 
which have temperature and density distributions similar 
to thos e of AFGL 2591. 

The lPotv et al.l l(2002l) model has most gas-phase sul- 
phur initially in S at T < 100 K (S/H2=6x IQ-^) and 
in H2S at T > 100 K (H2S/H2=1.6x IQ-^) to mimic the 
effect of ice evaporation. This choice of initial conditions 
is consistent with the nondetection of the [S I] line at 
25.249 ^m toward our sources with ISO. The upper limit 
of T <G.01 in absorption implies A^(S) < 5 x 10^^ cm~^, 
so that for 7V(H2) - 10^3 cm"^ S/H2< 5 x IQ-^ in- 
dicating that not all gas-phase sulphur in the inner en- 
velope can be in atomic form. The inferred abundances 
imply that SO2 is one of the dominant sulphur-bearing 
gas-phase molecules in the warm inner envelope. However, 
the sum of all sulphur-bearing molecules in the gas phase 
is still much lower than the value derived for diffuse clouds 
(Savage & Sembach 1996). Unless a major sulphur-bearing 
gas-phase species has been missed in this survey, ~90% of 
the sulphur must be in a solid form which has not yet been 
detected. 

Infrared observations of our sources with ISO-SWS do 
not support the assumption that II2S is the main sulphur 
reservoir of the grain mantles. Toward W33A, a source 
with large ice abundances, the 3.98 /zm band of solid II2S 
is not detected to r < 0.02, which using a band strength 
of 5 X 10-18 cm mol-i gives iVs(H2S) < 3.6 x 10^^ cm'^, 
or A/'s(H2S)/iVs(H20) < 0.03 (W. Schutte, priv. comm.). 
However, solid OCS was detected with an abundance 
of 7V.(0CS)/A^.(H, 0) = (0.4 - 1) x IQ- ^ or OCS/H2 
= (4 - 14) X 10"*^ l|Pahimbo et al.lll997^ . Evaporation 
of solid OCS may explain our measured gas-phase OCS 
abundances which are otherwise unaccounted for. W 33 A 
is the only source in our sample for which both gas-phase 
and solid OCS have been detected. The results imply 
a gas/solid ratio of ~0.5, much higher than that found 
for e . g. H2O and CO2 in this source ijBoonman et alJ 
l2003t Ib oonman fc van Dishoeckll2003|) . Observations of 
both gas-phase and solid-state OCS toward other sources 
are needed to determine the role of grain-mantle evapora- 
tion. Observations of OCS lines from a wide range of en- 
ergy levels may reveal an abunda nce increase at 1 00 K, 
as indeed found for IRAS 16293 l|Schoier et al.ll2002() . The 
low Tox of H2S and the high value of OCS ('!) I3.2|I support 
these conclusions. 

The effect of evaporating OCS-rich ice mantl es on hot 
core chemistry is explored bv IPotv et al.l l(2003t) . Briefly, 
the gas-phase abundances of CS, HCS"*" and OCS increase 
about linearly with the fraction of solid OCS, while those 
of H2S, SO and SO2 decrease about linearly. The H2CS 
molecule is not affected much. Since the gas-phase chem- 
istry of CS, HCS+ and OCS depends only weakly on time, 



the use of sulphur molecules as chemical clock diminishes 
when H2S is only a minor sulphur carrier in the ice man- 
tles. However, the 'jump' models for SO2 (Table indi- 
cate that either H2S or SO2 is present in the grain mantles. 
The limits on s olid SO2 from infrare d observations are not 
very stringent l|Boogert et alJll998(l . 

We ha ve compared our observations to the models of 
iDotv et~\. (2002), modified to have OCS in the evapo- 
rating ice mantles (S. Doty, priv. comm.). This model re- 
produces our observed abundances of CS, SO and H2CS 
in the outer envelopes (Table IHJl to within factors of a 
few for chemical ages of ~10^ yr. The abundances of SO2 
in the inner envelopes (Table are also reproduced for 
the same chemical age. The model abundances of HCS+ 
and SO2 in the outer envelopes are factors of ~10 below 
the observed values, which may be due to the adopted ini- 
tial conditions, or the chemical net work, which is based o n 
the UMIST reaction rate database l)Le Teuff et al.l2000() '^. 
The outer envelope abundances of H2S and OCS are un- 
derproduced even more (factors of ^100). However, the 
small number of observed lines forced us to assume con- 
stant abundances for these species, while they may in fact 
be confined to the inner envelopes. 

5.3. The role of shock chemistry 

As an alternative to ice evaporation models, we consider 
the effect of non-dissociative shocks on the molecular en- 
velopes surrounding the protostars. Such shocks occur by 
the interaction of the molecular outflows of protostars 
with their envelopes. Shocks fast enough to dissociate H2 
{v ^ 40 km s~i) would have a major impact on their molec- 
ular composition which our data do not indicate. Although 
our sources show outflows faster tha n that, both in in- 
frared absorpt ion ( M itchell et al .I1991I) and in submillime- 
ter e mission l)Mitchell fc Hasegawal Il99ll iMitchell et all 
I1992I) of CO, their flUing factor must be small. 

Our observed SO, SO2 and H2S abundances in the 
outer envelopes are consistent with those in postshock gas 
if sul phur was initially in molecular form ijLeen &: Grafj 
Il988i) . The observe d simila rity of t hese abundance s is not 
reproduced by the IPineaii des Forets et al.l ljl993l) model 
which has most sulphur initially in atomic form. However, 
the two models are hard to compare because one is a two- 
fluid treatment, while the other treats neutral and charged 
particles as one fluid. Before concluding that postshock 
gas reproduces our data, we recommend that a two-fluid 
calculation is performed with most sulphur initially in 
molecular form. 

Two other molecules suggest that the envelopes of 
young high-mass stars may have been processed by 
shocks, although neither one conclusively. First, ISO- 
SWS data show low ratios of gas-phase to solid-state 
CO2 r < 0.1: Ivan Dishoeckl 11991 and the CO2 gas phase 
abundance remains low through the 100-300 K temper- 
ature regime. After evaporating off grains, CO2 must 

^ http://www.rate99.co.uk 



12 



van der Tak et al.: Sulphur chemistry around massive young stars 



be promptly des troyed, which shocks can do in re- 
actions with H llCharnlev fc KaufmanI 2000) or pcr- 
haps Ho llDotv et al.ll2002HTalbi k Herbstll2002|) . Second, 
iHatchell fc Vitil tOO^ nieasured NS/CS ratios of 0.02- 

0. 05 in a sample of six hot cores and interpreted these as 
evidence for shocks. The main reactions to form NS require 
SH and NH which are produced from OH. Shocks use OH 
to form H2O and suppress the production of NS. The val- 
ues of NS/CS=0.001 - 0.01 measured here may indicate 
tha t shocks pla y a role too. However, for i = 3 x 10"* yr, 
the lDotv et aD j2002) model also predicts NS/CS - lO^^, 
so this ratio cannot be used to demonstrate the influence 
of shocks. 

6. Conclusions 

Submillimeter observations of SO, SO2, H2S, H2CS, OCS, 
NS and HCS"*" toward nine embedded massive stars show 
that: 

1. High-velocity gas contributes 10-50% to the CS, SO 
and SO2 emission in 15 — 20" beams. This gas may 
either trace the outflows which are also seen in CO, or 
they may trace infall motions. 

2. The SO2 abundance increases from dark cloud levels 
in the outer envelope (T < 100 K) to levels seen in hot 
cores and shocks in the inner envelope (T > 100 K). 

3. Molecular abundances are consistent with a model of 
ice evaporation in an envelope with gradients in tem- 
perature and density for a chemical age of ^10^ yr. 

4. The high observed abundance of OCS, the fact that 
Tex(0CS)»Tex(H2S), and the detection of solid OCS 
and limit on solid H2S all suggest that OCS is a major 
sulphur carrier in grain mantles, rather than H2S. 

5. For most other sulphur-bearing molecules, the source- 
to-source abundance variations by factors of up to 10 
do not correlate with previously established evolution- 
ary trends in temperature tracers. These species probe 
the chemically inactive outer envelope. Our data set 
does not constrain the abundances of H2S and SO in 
the inner envelope, which, together with SO2, are re- 
quired to use sulphur as a clock. 

These conclusions could be followed up by, respec- 
tively: 

1. interferometric observations at ^ 1" resolution of SO 
lines from a wide range of energy levels. 

2. similar interferometric o bservations of SO2. 

3. updating the model of iLeen fc Grafj l)l988j) using a 
two-fluid treatment, and extending it to later times. 

4. a multi-line study of OCS, covering energy levels well 
above and well below 100 K. 

5. observations of high-excitation SO and H2S lines, e.g. 
SO 430 and 473 GHz and H2S 369 GHz. 

Items 3 and 4 are within the capabilities of current 
computational and observational facilities. The CSO and 
JCMT may be able to carry out item 5 in very good 



weather, but the future Atacama Pathfinder Experiment 
(APEX) can benefit from an even better site. Current in- 
terferometers are not sensitive enough for items 1 and 2, 
which require observing several lines and several sources. 
However, future instruments such as the Combined Array 
for Research into Millimeter Astronomy (CARMA) and 
the Atacama Large Millimeter Array (ALMA) will be well- 
suited for these projects. 
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